Mean-field solution of the Potts glass near the transition temperature to the ordered 

phase 
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We expand asymptotically mean-field solutions of the p < 4 Potts glass with various levels of 
replica-symmetry breaking below the transition temperature to the glassy phase. We find that 
the ordered phase is degenerate and solutions with one hierarchy of spin replicas and with the full 
continuous replica-symmetry breaking coexist for p > p* ~ 2.82. The latter emerges immediately 
with the instability of the replica-symmetric one. Apart from these two solutions there exists also 
a succession of unstable states converging to the solution with the continuous replica-symmetry 
breaking that is marginally stable and has the highest free energy. 

PACS numbers: 64.60.Cn,75.50.Lk 



I. INTRODUCTION 

Models of spin glasses have been attracting a consid- 
erable attention of condensed matter theorists for more 
than three decades. Original motivation for construct- 
ing microscopic models of spin glasses came from ex- 
perimental observation of an unusual low-temperature 
behavior of magnetic impurities randomly diluted in a 
nonmagnetic mctali It soon became clear that a new 
type of long-range order of magnetic impurities emerges 
in spin glasses, qualitative understanding of which de- 
mands theoretical modeling. Rather soon Edwards and 
Anderson proposed a generic model of Ising spins with 
a randomly distributed exchange interaction^ A naive 
mean-field solution^ now known as the Sherrington- 
Kirkpatrick (SK) model, started a wave of intense inves- 
tigation of Ising and Hcisenbcrg spin glasses. The reason 
for the extended interest in spin-glass theories was in- 
sufficient understanding of the inconsistency of the SK 
mean-field approximation. A fully physically consistent 
solution was then proposed by Parisi via the replica trick 
used to handle quenched averaging over random configu- 
rations of the spin exchanged A rigorous proof of exact- 
ness of the Parisi free energy of the SK model, known as 
a full replica-symmetry breaking (FRSB), was completed 
only recentlyi££ 

The full mean-field solution of the Ising spin glass is 
not only complicated in its analytic structure but also 
in its physical interpretation^ That is why simpler mod- 
els, random energy^ p-spin^ Potts^ or quadrupolar— 
glass models, have been introduced so that the origin and 
the meaning of replica-symmetry breaking in frustrated 
spin models can be better understood. All four mod- 
els behave differently from the Hciscnberg spin glasses. 
They show temperature intervals where a first step to- 
ward the Parisi solution in the replica trick, the so-called 
one-level replica symmetry breaking (1RSB), is locally 
stable. While 1RSB is the true equilibrium state of the 
first model in low temperatures, it is locally stable only 
in a small interval of temperatures close to the transi- 
tion temperature to the glassy phase in the latter three. 



A full replica-symmetry breaking solution with infinite- 
many hierarchies of replicated spin variables should lead 
to a marginally stable equilibrium state at very low tem- 
peratures there. The genesis of FRSB and the way a 
solution with a continuous distribution of hierarchies of 
replicated spin variables is reached in these models with- 
out reflection symmetry in spin variables has not yet been 
fully clarified. 

We recently demonstrated that, although the first level 
of replica-symmetry breaking below the transition tem- 
perature to the glassy phase in the mean-field model of 
the Potts glass is locally stable, a full continuous replica- 
symmetry breaking solution coexist with it and has a 
higher free energy^ This result suggests that FRSB 
leads even in models without reflection symmetry to the 
true equilibrium state independently of the fact that a 
replica-symmetry breaking solution with a finite num- 
ber of replica hierarchies is locally stable. It has been 
assumed that the Parisi continuous replica-symmetry 
breaking in models without spin-reflection symmetry fails 
or is inconsistent ! 13 ' 14 A natural question arises when 
a locally stable discrete and marginally stable continu- 
ous RSB solutions coexist: How does the FRSB solution 
emerge when a discrete RSB solution with finite-many 
replica hierarchies no longer decays to solution with a 
higher number of replica hierarchies. The full RSB state 
has been assumed to emerge only below the temperature 
T-2 at which the 1RSB solution becomes again unstable. 
Heuristic arguments were used to suggest a cascade of 
successive transitions below the instability of 1RSB i 15 i 16 

The aim of this paper is to demonstrate explicitly the 
genesis of the Parisi solution with a continuous FRSB 
in the model of the Potts glass with p states. We use 
an asymptotic expansion of the Parisi solution near the 
critical temperature and show that for p < 4 the solution 
with FRSB emerges at the transition point at which the 
replica-symmetric solution gets unstable and the system 
undergoes a continuous transition to the glassy phase. 
We find that near the transition temperature to the 
glassy phase there exists apart from a locally stable 1RSB 
solution for p* < p < 4 also a series of unstable solutions 
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with K = 1,2,... discrete hierarchies breaking replica 
symmetry and converging towards a marginally stable 
Parisi solution with a continuous FRSB. 



II. MEAN-FIELD MODEL 



A. Potts Hamiltonian and a replicated mean-field 
solution 



The Potts model is a generalization of the Ising model 
to more than two spin components. The original formula- 
tion of Potts^ with Hamiltonian H p = — Yli<j Jij^ni.nj 
where n% = 0, ... ,p — 1 is an admissible value of the 
p-component model on the lattice site Ri, is unsuitable 
for practical calculations. The Potts Hamiltonian can, 
however, be represented via interacting spins* 8 



i^./,.S,.S, ^l.-H, 



(1) 



'■J 



where S, = {sj, . . . sf" 1 } are Potts vector variables tak- 
ing values from a set of state vectors {ex}x=i- Functions 
on vectors e A are in equilibrium fully defined through 
their scalar product 



V 



-1, 



(2a) 



where a e {1, ...,p— 1}. We use the Einstein summation 
convention for repeating Greek indices of the vector com- 
ponents indicating a scalar product of the Potts vectors. 

The components of the Potts vectors obey the following 
sum rules 



A=l 
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from which we can construct their explicit representation 







£a = 



p(p-«) 

p+1 — a 
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A < a 
A = a 

A > a 



(3) 



It is evident from the representation of Potts vectors in 
Eq. Q that the Potts variables are not symmetric around 
zero for p > 2, that is the Potts model does not possess 
spin-reflection symmetry. In case of p = 2 the Potts 
model reduces to the Ising model. 

Frustration due to a quenched randomness in the Potts 
model is introduced via static fluctuations of the inter- 
action parameters Jij being distributed randomly with 
probability 



P(Jij) = 



1 



(2ttJ 2 /N) 1 /- 



exp ■ 



-(Jij - JqY 
2.P/N 



(4) 



where Jo = iV _1 J^j Joj i s the averaged (ferromagnetic) 
interaction and N is the number of lattice sites. Random- 
ness is assumed to be quenched, that is, the equilibrium 
free energy is averaged over the configurations of the spin- 
spin interaction. The spin-spin interaction is long-range 
with an infinitesimal (-/V -1 ) variance of its fluctuations 
(mean-field model). 

The standard way to derive a mean-field approxima- 
tion to frustrated models with random exchange interac- 
tions is to use the replica trick with which a quenched 
averaging is transformed to an annealed one of a repli- 
cated system^ There is, however, a price we pay for this 
simplifying transformation. We must perform a limit to 
zero number of replicas, which demands analytic con- 
tinuation of the solution in the replica number. This is 
possible only under specific conditions with an appro- 
priate symmetry of the order parameters in the repli- 
cated phase spaced There is a possibility to avoid the 
replica trick with the limit to zero replicas in that we de- 
mand thermodynamic homogeneity of the averaged free 
energy^ Although replicas are used, there is no need for 
the limit to zero number of replicas. Instead, indepen- 
dence of the replication index is demanded in the real- 
replica approach. The replicated mean-field free energy 
of the Potts glass with the gaussian distribution of spin- 
spin interactions and the probability distribution from 
Eq. ||1J) is 



U = 



/3J 2 1 
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where v is the number of replicas, h a = h a + Jom" + 
Jy a y/q is an effective magnetic field, and rj a are Gaussian 
random fields with a {p — l)-dimensional measure 



p-i 



V 



(p-i) 



(y) = n 



drf 



exp 



(y a Y 



We further denoted an effective "ferromagnetic" ex- 
change Jo = Jo + PJ 2 (p — 2)/2. The order parameters 
in the replicated free energy f v are the averaged square 
local magnetization q = N^ 1 mfmf and local over- 



lap susceptibilities x"b = ^ 1 S 
for a 7^ b measuring the linear response of the replicated 
system to a small inter-replica interaction. If the free 
energy is thermodynamically homogeneous, the overlap 
susceptibilities must vanish in equilibrium, a saddle point 
of the replicated free energy. Free energy is then inde- 
pendent of the replica index v. It reduces in this case 
to the replica-symmetric one having a representation for 



s?„ s. 
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the isotropic Potts model (zero magnetic field and no 
long-range ferromagnetic order, mf = 0) 

-Mrs = ^( P -l)(q-l) 2 + J 2Vi)(y)lnZ (y) • ( 6a ) 



We denoted the local partition function of the Potts 
model in a random magnetic field 



z (y) = E exp{/V9y a el} 



(6b) 



and {\}i 
reads 



, Xk ■ The initial zero- level partial sum 



Z (y,X) = Z^(y,{\} 1 ) 



A=l 



We need to break the replica symmetry so that to check 
thermodynamic homogeneity of the equilibrium free en- 
ergy from Eq. ©. That is, to test stability of the replica- 
symmetric solution with respect to replications of the 
phase space of the relevant order parameters. 



B. Free energy with discrete hierarchies of 
replica-symmetry breaking 

A natural way to start with replications of the original 
model is to use simply two replicas. It was done in detail 
for the Ising spin glasses It was demonstrated there that 
although free energy was lowered in the system with two 
replicas the instability and thermodynamic inhomogene- 
ity of the replica-symmetric solution was made worse. It 
is necessary to continue analytically the replicated free 
energy to an arbitrary positive replication index to ana- 
lyze dependence of free energy on the replication index. 
Making the replication index a continuous variable is pos- 
sible only if the symmetry of the matrix of the local over- 
lap susceptibilities possesses the symmetry introduced by 
Parisi in his construction within the replica trick. Instead 
of numbers of replicas used we introduce a number of hi- 
erarchies of replicas that distinguish different solutions. 
A solution with K hierarchies is determined from a sad- 
dle point of a free-energy functional 



Pf K (q,{A X i},{mi}) 



3 2 J 2 



(P-1) 1-5 



K 



K 



1=0+1 



y© (p _ 1) (y)ln^(y) (7a) 



The equilibrium state for this free energy is character- 
ized by the averaged square magnetization q and a set of 
K pairs Axi, mi for I = 1, 2, . . . , K. Thermodynamic ho- 
mogeneity of a isT-level hierarchical solution is achieved if 
Axk+i = 0, which leads to independence of the system 
on the next replicating parameter mK+i^ The overlap 
susceptibilities < A\i < 1 form generally a decreasing 
sequence, since we demand that the last one should van- 
ish in the thcrmodynamically homogeneous system. The 
indices counting the replica hierarchies m; can be arbi- 
trary. There is, however, a degeneracy in the hierarchical 
free energy. We obtain fx+i = /if if mK+i = 0,mfe,oo. 
It means that we obtain at least one new equilibrium 
solution for a pair Axrt+i, mK+i with mx+i < rax- 
The new mean-field solution is acceptable if it leads to 
a smaller thermodynamic inhomogencity measured by 
Axk+i < A%if. In the asymptotic solution near the 
critical temperature of the Ising spin glass the new so- 
lution with rriK+i < tuk leads to a higher free energy 
Jk+i > /if 21 ' 22 There is also another stationary solu- 
tion for rriK+i > ttik that, in the Ising model, lowers the 
free energy and worsens thermodynamic inhomogencity. 
Hence, it is unacceptable. We can do a similar analy- 
sis of the XRSB free energy, Eq. (J7J), near the transition 
temperature of the Potts glass. 

We start with a solution with a first level of replica- 
symmetry breaking. An explicit mean-field free energy 
of the Potts glass with 1RSB reads 



/?/i = ^(p-l) 



(1 - q - A X f - mA X (A X + 2q) 



+ - 1 2Vi)(y) ln / v (p-dW 4 P) (y.A) 



(8) 



where 



Zf(y,Wi+i) = 



V 



(AO(^) (y,{A},) r 



The partition function for the p-state Potts model is con- 
(7b) structed by using the representation in Eq. ©. For the 
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three-state model we obtain explicitly 



K 



4 3) (y. Wi) = exp { /Wv^ U^y 1 + £ V^O^' 



3=1 



exp < (3 J 



exp < — j3J 



(9) 



Properties of the mean-field theory of the Potts glass 
with 1RSB have been analyzed by several groups i 15 ' 23 i 24 
The most prominent feature of the solution of the Potts 
mean-field model with 1RSB and the number of states 
p > 3 is its local stability near the transition to the low- 
temperature glassy phase. One could conclude from this 
result that the hierarchical construction of the mean-field 
free energy stops just at 1RSB and no Parisi solution with 
FRSB exists in the region of local stability of 1RSB. Wc 
demonstrate on the asymptotic solution below the transi- 
tion temperature to the glassy phase that the mean-field 
equations of the Potts glass are degenerate and allow for 
a cascade of coexisting metastablc sates including the 
Parisi FRSB solution that is marginally stable alike in 
the Shcrrington-Kirpatrick model. 



C. Parisi solution with a continuous 
replica-symmetry breaking 

It is not necessary to derive a mean-field free en- 
ergy with FRSB via checking stability of free energies 
with finite-many hierarchies of discrete replica symmetry 
breakings. It is sufficient to look at the behavior of the 
hierarchical free energy in the continuous limit. It obeys 
a differential equation derived first by Parisi for the Ising 
spin glassi 4 - If we introduce a parameter A € (0, 1) and de- 
note X = qsA — q = J2i A-Xh then the A-dependent freee 
energy g(X, h) in the continuous limit of the replica sym- 
metry breaking hierarchy with K — > oo and Axi oc K^ 1 
obeys a differential equation that for a many-component 
spin model reads 



dg(X,h) 
dX 



X 
~2 



d 2 ff(A,h) 
dh a dh a 



m(A) 



dg(X,h) dg(X,h) 



dh° 



dh c 



(10) 



where m(A) is a continuous limit of the replica indices 
mi from the solutions with discrete replica hierarchies. 



Having this differential equation we can try to resolve it 
on a phase space of the order-parameter functions m(A). 
One of us recently suggested an explicit representation 
for the free energy obeying the Parisi differential equation 
(fT0|) 3£ It can easily be generalized also to the Potts glass 
for which we obtain 







/•i 



{p-l)(l- q -Xf 



2 [p-l)X J dXm{X)[q+X(l-X)] + (g(l,y^q)} y 



(11a) 



where (F(y)) y = / P (p _ 1} (y)F(y) and 
' X 



g(v, h) = T A exp 



2 Jo 

+m(X)g' a (X,h + h)d h 



dX [d ha d ha 

.go(h + h)| fi=0 (lib) 



with the initial local free energy #o(h) = 
ln^A=i cxp{f3h a e^}. Wc introduced an evolution 
operator represented via a "time-ordering" operator 
T\ ordering products of A-dependent non-commuting 
operators from left to right in A-dccreasing succession. 
We further denoted g' a (X,h) = dh a g{X,h). We used an 
auxiliary vector field h = (h 1 , h 2 , . . . hP^ 1 ) to generate 
the necessary derivatives of the bare free energy go. 

Free energy, Eq. (fTTj) . is stationary with respect to vari- 
ations of variables q and X and function m(A) for every 
A € (0,1). It is straightforward to show that g(X,h) 
obeys differential equation (fTO]) . Although representa- 
tion (jlip of a mean-field free energy with a continuous 
FRSB of the Potts glass is implicit, it is self-contained 
and allows us to derive explicit stationarity equations 
for the order parameters and other quantities. It also 
enables to reach directly approximate or asymptotic so- 
lutions without using solutions with finite numbers of 
replica hierarchies and studying their stability. We use 
representation (fTT]) here to derive an asymptotic solution 
with a continuous replica-symmetry breaking of the Potts 
glass with p < 4 states below the transition temperature 
to the glassy phase. 



III. ASYMPTOTIC SOLUTION NEAR THE 
CRITICAL TEMPERATURE 

Full-scale solutions of the mean-field equations of the 
Potts glass are not available due to their complex struc- 
ture. What we can explicitly obtain are only asymptotic 
limits of the solutions with various degrees of replica- 
symmetry breaking. We hence expand asymptotically 
the replica-symmetric solution, the solutions with dis- 
crete finite-many hierarchies of the replicated spin vari- 
ables as well as the Parisi solution with the continuous 
order-parameter function. The expansion coefficients will 
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be calculated with program MATHEMATICA. We ex- 
pand the corresponding free energies to the fifth order 
in the small expansion parameter being r = (1 — T/T c ), 
where T c = J = 1 is the temperature of a continuous 
transition to the low-temperature glassy phase at which 
different mean-solutions can be distinguished. Since we 
assume a continuous transition to the glassy phase, our 
analysis restricts to the Potts model with the number of 
states p < 4. The derived asymptotic solutions allow us, 
however, to analyze the asymptotic behavior of the glassy 
phase below the transition temperature as a function of 
a continuous parameter p. 



A. Replica-symmetric solution 

The simplest mean-field state is a stationary solution of 
the replica-symmetric free energy from Eq. ^ . The only 
order parameter is the average of the square of local mag- 



crgy 



nctizations q — N 1 £\ 



This pa- 



rameter vanishes in the high-temperature paramagnetic 
phase and starts to grow continuously from zero below 
the transition point to the glassy phase. The correspond- 
ing stationarity equation for the replica-symmetric order 
parameter derived from free energy, Eq. ([6|), reads 



(p-l)q+l =p 



p f 

•E / I3 (p-i)(y) 

A=1 J 



exp{2f3^qy a e^} 

Zo(y) 2 



(12) 



We do not want to evaluate fully the gaussian integrals on 
the right-hand side of Eq. (| 12[) but rather only near the 
transition temperature. Since the transition is expected 
to be continuous we can assume that the order parameter 
q is small and expand free energy, Eq. ([5]), into a power 
series in q. We cut the expansion at the fifth order. The 
explicit expression for the expanded free energy is given 
in Appendix in Eq. (|Alj) . The first term in the expansion 
of the free energy proportional to q 2 changes sign at /3 = 1 
indicating a continuous transition to an ordered phase. 
The transition is continuous where the third order of the 
expansion is positive, that is, for p < 6. 

It is sufficient to expand order parameter q to the third 
order of the the expansion parameter r = 1 — T so that 
to obtain the expansion of free energy to the fifth order. 
The expansion calculated with program MATHEMAT- 
ICA reads 



Q = 



At 2(-7p 2 -60p+180)r 2 



6 — p 3(6 — p) 3 

8 (p 4 + 300p 3 - 1044p 2 + 4320p - 7776) t 3 
9(6 - pf 



(13) 



Inserting this expansion into Eq. (|A1|) we obtain an 
asymptotic expression for the replica-symmetric free en- 



8r 3 



Hp 2 



252 t 4 



RS 



p-l'"" 3(6 -p) 2 3(6 -p) 4 

2 (29p 4 + 1320p 3 + 54360p 2 - 294624p + 421200) r 5 

+ 45(6 -pf ' 

(14) 

The asymptotic expansion based on smallncss of the or- 
der parameter breaks down at p = 6 indicating a change 
in the way the transition to the glassy phase occurs. No- 
tice that the expanded free energy of the Potts glass in 
Eq. (fT4|) coincides with that of the Ising spin glass for 
p = 2& 



B. 1RSB solution 

The next step beyond the replica-symmetric solution is 
a state with the first level of replica-symmetry breaking. 
Its free energy is described by three order parameters 
q, Ax, m. We derive their defining equations from the 
stationarity point of the 1RSB free energy in Eq. (JS}. To 
be able to write down these equations in a condensed way 
we first introduce a useful notation. 

Z (y,A) = Y,E A (y,X) , (15a) 

A 

E A (X,y) = exp{/3 [^y a + \/A^A Q ] e a A } . (15b) 



Further on we will need 



t a = 



z (y,X) 
Z (y,X) m 
(Z (y,A)™) A 



(15c) 
(15d) 



An equation for the equilibrium order parameter q de- 
rived from dfi/dq = reads 



(p - l)(q + (1 - m)A X ) = (1 - m)((p t a t a ) x ) y 
+ m((p t a ) x (p t a ) x ) y . 

Analogously from dfi/dAx = we obtain 

(p-l)(q + A X ) = ((p t a e) x ) y . 

Finally the equation for parameter m is 

m^(p - l)A X (2g + A X ) = ((pin Z (y, A)) A >„ 

-(In^i(y)), 
[/2? (p _ 1) (A)Z (y,A)"f /m 



(16) 
(17) 



(18) 



where Zi (y) = 

(Zoix^r) 1 /- 1 . 

We simplify the stationarity equations in that we as- 
sume that parameters q and Ax are small near the tran- 
sition to the glassy phase. We then expand free energy 
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/i from Eq. ([8]) to the overall fifth order in these parame- 
ters. An explicit expression for the expanded free energy 
is given in Appendix, Eq. (|A2[) . It allows us to derive ex- 
plicit equations for the order parameters of the 1RSB so- 
lution. They are straightforward to derive but are rather 
lengthy. That is why we do not list them explicitly. We 
give here only the result of the asymptotic expansion in 
the small expansion parameter r. The results were again 
derived with the aid of program MATHEMATICA. The 
glassy phase of the Potts model allows for multiple solu- 
tions with discrete replica-symmetry breakings. 

We found a double degeneracy of the 1RSB solution. 
Both solutions have the same parameter m to the two 
lowest orders in r 



. p - 2 36 - 12p + p A 

m = 1 ; ; T 



8(4 -p) 



(19) 



One non-trivial 1RSB solution then leads to order pa- 
rameters 



'1 



«=0, 



(i) 



228 - 96p + p 2 2 



A-p 6(4 -p) 3 

while the second one has both parameters nonzero 



(20a) 
(20b) 



(2) ^ -12 + 24p-7p 2 2 



A X (2) = 



3(4-p)*(p-2) 
2 360 - 204p - 6p 2 + 13p 3 2 



6(4 -p) 3 (p- 2) 



(21a) 
(21b) 



Both the solution have the same asymptotic free energy 
to the fifth asymptotic order 



p-1 



\IRSB 



3(4 -p) 



(p(llp - 102) + 204)r 4 (p(p((18744 - 1103p)p - 120648) + 325728) - 317232)r 5 



12(4 -pf 



720(4 -p) 5 



(22) 



Unlike the replica-symmetric solution we can see that the 
asymptotic expansion with small parameters q and Ax 
breaks down already at p = 4 above which we expect 
a discontinuous transition from the paramagnetic to a 
1RSB state at T > 1. 

The 1RSB solution has a higher free energy than the 
replica-symmetric one. The difference is of order r 3 , 



1RSB 



fas = 



(p-2) 2 (p-l)r 
3(4-p)(6-p) 2 



(23) 



The two stationary states of the 1RSB free energy, 
Eq. (|A2p . behave differently clS iX function of the pa- 
rameter p. The former solution is physical for all val- 
ues of p unlike the latter that becomes unphysical for 
p > p* w 2.82 where q^ 2 ' from Eq. ()21a|) turns negative. 
It is also the region of the parameter p where the first 
solution is locally stable as can be seen from the stability 
function 



A< 0) = p - 1 - /3 2 ( (p [t aP " (1 - m)t a t a ] 



XI y 



r 2 (p-l) 



6(4 -p) 2 

that is positive in this region. We denoted 

Zo(y,A) 



(7p 2 -24p+12) . (24) 



(25) 



The first 1RSB solution q( x \&x is not physically 
inconsistent and is locally stable for p > p* . It should not 
break into solutions with a higher number of hierarchies 
of replicated spins. That is why it has been considered 



I 

as the equilibrium state of the Potts glass on a "high- 
temperature" interval of the glassy phase. 

The two asymptotic solutions behave differently when 
approaching the Sherrington-Kirkpatrick model, p — > 2. 
The first is regular and goes over at p = 2 to a new 1RSB 
solution with q = 0. Since both A% and m are nontrivial, 
it is a solution breaking replica symmetry. This solution 
was not discussed in Ref. Wi as we excluded q = case. 
The second solution is singular in the limit p — > 2. The 
singularity can be seen first the second order in the ex- 
pansion parameter r. Consequently, the asymptotic lim- 
its t — > and (p— 2) — > do not commute and the result 
depends on the value of ratio rj (p — 2). Due to this sin- 
gularity we cannot unambiguously continue the second 
solution the Sherrington-Kirkpatrick model. This non- 
analyticity is connected with emergence of spin-reflection 
symmetry in the two-component Potts (Ising) model. 



C. A'RSB solution 

Although we have a locally stable solution of the 1RSB 
free energy for p > p* , there is another unstable solution 
that decays to a solution with a higher number of hier- 
archies of the replicated spins. We hence can investigate 
possible solutions of free energies with an arbitrary num- 
ber of spin hierarchies. 

We start with stationarity equations for the i^RSB 
free energy from Eq. (JT)). The equation for the order 
parameter q reads 

(p-l)q = ((t a ) K {t a ) K ) y , (26) 

where we denoted (F)i(y, {A}; + i) = 

/2Vu(A,) Pl ( y ,{^}i)(F)i-i(y,{Wi), Pi{y,{Mi) = 
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Z i _ 1 (y,{A} i ) m V(^ I -i)A ! (y,{A}i+i). Wc further abbre- 
viated (X)\ t = J 2?( p _!)A; X. Equations for the overlap 
susceptibilities are 

(p-l)A XJ = (((t«)j-i(t )i-i>Jc)» 

- (((t a )l(t a )l}K}y ■ (27) 

They are accompanied by equations for scaling parame- 
ters mi having the following form 

with I = 1, ...,K. 

One cannot solve these equations fully but we can 
reach their solution in the asymptotic limit below the 
transition temperature to the glassy phase as done for 
the replica-symmetric and 1RSB solutions. We obtain 
only a single solution for K > 1 within the leading-order 
asymptotic limit r — > 



q 

Ax 



K 



1 12-24p + 7p 2 2 
~3K 2 {4- P ) 2 (p-2) T 
1 2 



,K 



K{A-p) ' 
p- 2 2 



4 — p 



3 + |p 



7 r 

3 — 6p + —p' 



P 

21-1 



2K 



(29) 
(30) 

(31) 



We can see that the ifRSB solution behaves unphysi- 
cally in the same way as the second 1RSB solution does. 
The averaged square of local magnetizations is negative 
for p > p* where the first 1RSB solution is locally sta- 
ble. Negativity of q means that local magnetizations are 
imaginary and the solution is unphysical. This deficiency, 
however, decreases with the increasing number of spin hi- 
erarchies and disappears in the limit K — > oo. It is anal- 
ogous to negativity of entropy in the low-temperature 
solutions of 7f RSB approximations of the Sherrington- 
Kirkpatrick model. 

Free-energy density of the ifRSB solution increases 
with increasing the number of hierarchical levels 



p-1 



KRSB 



3(p - 4) 



(p(Up- 102) + 204)r 4 
12(p - 4) 3 



{p{p{p{16p - 265) + 1686) - 4532) + 4408)i 
10(p-4) 5 

(7p 2 - 24p+ 12) 2 r 5 
720A" 4 (p - 4) 5 



(32) 



and reaches its maximum at K = oo. Since the overlap 
susceptibilities are linear in 1/K, the limit K — > oo leads 
to a solution with a continuous order-parameter function. 



D. Continuous replica-symmetry breaking 

We showed in the preceding section that the Parisi 
continuous full RSB solution is not isolated and there is 
a cascade of discrete -ftTRSB states, even in the region of 
stability of a 1RSB solution, converging to the full contin- 
uous RSB. The asymptotic expansions of these solutions 
are singular in the limit of the Ising spin glass, p = 2. The 
singularity, however, vanishes in the limit K = oo and the 
asymptotic expansion in r of the Parisi solution appears 
to be analytic around p = 2. The leading asymptotic 
order of the full RSB solution can then be obtained from 
the limit of the discrete ifRSB approximations of the pre- 
ceding section. To lift the degeneracy in the asymptotic 
region T — > T c , we have to expand the order-parameter 
function to higher powers of the small parameter r. We 
use the explicit representation in Eqs. (fTTj) to determine 
the asymptotic limit of the Parisi solution near the tran- 
sition temperature. 

The characteristic function for the Parisi solution is a 
"dynamical magnetization" with the following represen- 
tation via an evolution operator 

^(A,h)=E(X,h;A,0)o 5 ' 



T„ exp {X J dv 



-d ni3 d hf , +m(%a(i',h + h)<9 s 



I," 



%)(h + h) 



dh a 



(33) 



h=0 



The defining equation for the stationary value of pa- 
rameter q reads 



P 2 (p-l)q=(g' a (l,h ri )g' a 0-,h n )), 



(34) 



where = h + rj^/q. Parameter X is determined from 

(3 2 (p - 1)X = (E(X,h v ; 1,0) o [g'o.Mg^M]), 

- (g' a (l,h v )g'M,S))v (35) 

Order-parameter function m(A) is obtained from an iden- 
tity 

/3 2 (p-l)XA = (E(X,h ?) ;l,0)o [g'oJh^M]^ 
- (E(X, h„; 1, A) o [g' a (X, h v )g' a (X, h,)]), (36) 

valid for A G (0,1). 

We next use the knowledge from the limit K — > oo of 
the discrete if RSB solutions, namely 



<?;(!, 0) =0 



(37) 



expressing consistency of the continuous limit. 

Equation (|36|) holds for all index variables A E 
(0, 1). Its form is unsuitable for determination of order- 
parameter function m(A). We can, however, perform an- 
alytic operations on both sides of this equation so that to 
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transform it to a more suitable form. Applying deriva- 
tive with respect to A leads to a condition of marginal 
stability.— It reads 

f3 2 (p -l) = (E(X,0;l,X)o [g'^(XM p (X,0)] . (38) 

To derive this form we used representation and dif- 
ferential equation (fTO]) . Further derivative of the above 
equation with respect to A leads to an explicit represen- 
tation for the order-parameter function 



E(X,0;1,X) 


o 






E(X,0;l,A)o 





Evaluating this expression at the transition point T c = 1 
where X = and E = 1 we obtain 2m(0) = p — 2. Below 
the critical temperature where X > we expand order- 
parameter function m(A) in XX and X as independent 
parameters 

m(A)= m\j,k]X j X k (40) 

j=0,k=j 



where k is the order of the asymptotic expansion and 
determines the asymptotic precision. We expand anal- 
ogously the A-dcpcndent free energy so that to be able 
to resolve evolution operator E. We have to keep depen- 
dence of free energy on an external magnetic field and 
hence 

5 (A;h)= 50 (h)+ J2 9\j,k;h}X j X k . (41) 

j=i,k=j 



We expand all quantities in powers of X and AX and 
solve each equation for individual orders independently. 
We do not list here all equations for the expansion pa- 
rameters. Nontrivial expansion parameters m[j, k] and 
9(j> k) = flL?) 0] to the order k = 5 are listed in Ap- 
pendix. 

The low-temperature glassy phase is reached when 
there is a nontrivial solution of Eq. f|35[) for parameter 
X. In the asymptotic limit X — > we obtain an equa- 
tion allowing for a nontrivial solution if /3 > 1 



= ^X 3 /3 8 (/3 2 (4p 2 (3m[0, 0] - 19) + 6p (m[0, 0] 2 - 29to[0, 0] + 114) 

-12 (2m[0, 0] 2 - 33m[0, 0] + 105) +p 3 ) + 2(p - 4)(3m[0, 1] + m[l, 1])) 

+ ^X 2 f3 s (2p(m[0, 0] - 10) - 8m[0, 0] + p 2 + 50) + (p - 4) A/3 6 + @ 2 (f3 2 - 1) . (42) 

Using expansion coefficients m[j,k] from Appendix we obtain the first four exact powers in parameter t = (/? — l)//3 

2 (p(p+ 12) - 36) 2 (560- 4p(p(llp- 57) + 136)) 3 
-A- = T ; T ; tt t 

A-p (4-p) 3 (4-p) 5 

r 4 + 0(r 5 ) (43) 



(p{p(p{p(71p - 705) + 5832) - 25832) + 54960) - 44560) , 



(4-pY 

being small (finite) for p < 4. 

We can analogously evaluate free energy of the Potts glass. Using again expansion coefficients m[j, k] we obtain an 
asymptotic expansion up to the fifth order in X of g(l, 0) 

g(l,0) _ log(p) 
p — 1 p — 1 

+ ^f- + \(P - 4)A 2 /? 4 - -L(p(p + 38) - 112)X 3 /3 6 - g^(p(p(19p + 194) - 3828) + 7896)X 4 /3 8 

A 5 /3 4 (4m[3, 3] + 3(p(p((1765- 576p)p + 45394)- 322428) + 493032)/? 6 ) 6 
+ ^ + 0(X). (44) 



We insert the asymptotic values of parameters X and parameter r 
to [3, 3] and obtain an explicit dependence on the small 



-Mr) = + 



(p(llp - 102) + 204) 



p-V y ' 3(4 -p) 12(4-p) 3 
(p(p(p(16p - 265) + 1686) - 4532) + 4408) 5 

10(4 -p)5 T ' (45) 
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It is easy to demonstrate that free energy of the full RSB 
solution is higher than free energies of the discrete RSB 
solutions. We have 

P(Jc - Jkrsb) = - on „,,, (46) 



and 



P(f c - fus) 



720A' 4 (4 - p) 



(p-l)(p-2) 2 r 3 
3(4-p)(6-p) 2 



(47) 



Parisi-likc solution with a continuous order-parameter 
function has the highest free energy as in the SK model 
and represents the true thermodynamic equilibrium of 
the Potts glass for 2 < p < 4. 



IV. DISCUSSION AND CONCLUSIONS 

We studied in this paper the asymptotic behavior of 
the mean-field p-state Potts glass below the transition 
temperature to the glassy phase defined by instability 
of the replica-symmetric solution. We analyzed sep- 
arately the replica-symmetric solution, solutions with 
1RSB, A'RSB for K > 1, and a solution with contin- 
uous RSB. We separated the 1RSB scheme, since its 
free energy is the only one with two distinct station- 
ary states. We found that the RSB solutions peel off 
from the replica-symmetric one continuously for p < 4 so 
that the small expansion parameters are q, the averaged 
square of the local magnetization, and the overlap sus- 
ceptibilities A\i between neighboring hierarchies of spin 
variables I — 1 and I. We expanded these order param- 
eters to the fifth order in r = 1 — T/T c to distinguish 
individual states. The asymptotic expansion allowed us 
to analyze the behavior of the Potts glass also as a func- 
tion of parameter p as a continuous variable connecting 
solutions for 2 < p < 4. We were able to distinguish two 
regions where the solutions with finite-many replica hier- 
archies behave differently. The two regions are separated 
by a critical value p — p* « 2.82. 

The Potts model in the region with 2 < p < p* has 
only a single solution representing the true equilibrium. 
It is the solution with continuous RSB being marginally 
stable. The solutions for the schemes with finite-many 
replica hierarchies are all unstable as in the case of the 
Ising spin glass. The A'RSB solutions are, however, non- 
analytic around the Ising limit, p = 2. A singularity in 
the asymptotic expansion of these solutions emerges in 
the second order of r for p = 2. It means that there is 
no analytic continuation of the ARSB solutions from the 
Ising spin glass to the Potts one. This non-analyticity is 
caused by spin-reflection symmetry present only in the 
Ising model. The 1RSB free energy of the Potts model 
deserves special attention. Apart from the asymptotic 
solution with q =/= singular at p — 2 we found another 
solution with q = free of any singularity in the limit 
p — > 2. This solution can hence be analytically continued 



to p = 2 to a 1RSB solution with non-zero parameters 
Ax and m. Such analytic continuation has a lower free 
energy than the standard 1RSB solution with q > dis- 
cussed in Ref. [U It is interesting to note that the 1RSB 
solution of the SK model with q = decays to solutions 
with a higher number of replica hierarchies. The asymp- 
totic expansion of parameters A\ and m for the general 
A'RSB solution reads 



f 



(2(K 



(6A' : 



1) 



,K 



6A 3 
j + 1) 



1) 



o K 

(\2K 2 + 1)(2(A' — j + 1) — 1) 2 
6A' 3 

All these A'RSB solutions with q = are unstable, have 
a lower free energy than the A'RSB solutions with q > 
and converge towards the Parisi full RSB solution where 
q = as well. 

The RSB schemes with finite-many hierarchies change 
their behavior after passing a critical value p = p* where 
the stability function of the first 1RSB solution, be- 
comes positive, indicating its local stability. Moreover, 
the second 1RSB solution and all other ARSB schemes 
turn unphysical, since the averaged square of the local 
magnetization goes through zero and gets negative, q < 
for p > p* . It is, however, important that negativity of 
parameter q decreases with increasing A and approaches 
zero in the continuous limit A — > oo. The unphysical 
ARSB states converge then towards a solution with full 
continuous RSB for which q = 0. Hence, the contin- 
uous RSB solution does not experience any unphysical 
behavior for p > p* . It remains marginally stable for all 
values of p, has the highest free energy from all studied 
construction schemes, and is thermodynamically homo- 
geneous. The only observable change in the full RSB 
solution for p > p* , as discussed in Ref. [H, is a change 
in the sign of the derivative of order-parameter function 
m(A), cf. coefficient m[l, 1] in Appendix. 

To conclude, we analyzed the asymptotic behavior of 
the mean-field random Potts glass with the number of 
states p < 4 below the transition to the glassy phase. 
We demonstrated that a Parisi-like solution with contin- 
uous replica-symmetry breaking emerges simultaneously 
with the instability of the replica-symmetric solution. Its 
existence is independent of local stability of the 1RSB 
solution observed for p > p* rj 2.82. We found that the 
solution with continuous RSB in the Potts model is a 
limit of other unstable solutions with discrete RSB that 
decay towards it alike in the Shcrrington-Kirkpatrick 
model. These solutions arc unphysical, since the aver- 
aged squared local magnetization is negative, but the 
resulting limit with continuous RSB is free of any un- 
physical behavior. We studied in this paper only the 
Potts model with p < 4, but our construction of the 
asymptotic expansion in parameter X = qEA — qsK 
can be extended also to p = 4 + e for which I « 1. 
There we expect that replica-symmetry breaking solu- 
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tions emerge above the critical temperature of instability 
of the replica-symmetric one. The solution with contin- 
uous RSB does not seem to display any singularity at 
p = 4. The question to be answered in the Potts model 
with p > 4 is: which solution represents the true equi- 
librium state in the region of coexistence of the replica- 
symmetric and replica-symmetry breaking solutions? 

Research on this problem was carried out within 
project AV0Z10100520 of the Academy of Sciences of the 
Czech Republic. 



Appendix A: Expanded free energies near the 
critical transition temperature to the glassy phase 

Free energy of the replica-symmetric solution has only 
a single parameter q that we can use as the expansion pa- 
rameter below the critical transition temperature to the 
glassy phase. We expand free energy to the hfth order, 
although third order would be sufficient to distinguish it 
from the replica-symmetry breaking ones. We obtain for 
q < 1 



P 



jfRS = ^q 2 (/3 2 - 1) + ^ (P - 6 )l 3 + ig - 3 °P + 90 ) 1 + iio ( P - 1Up2 + U36p ~ 252 °) ? ■ ( A1 ) 



We explicitly expand only 1RSB from all discrete if RSB solutions where we use two small expansion parameters q 
and A%. It is necessary to expand it to the fifth order so that to distinguish it from higher-order RSB solutions and 
the continuous full RSB one. We obtain 

-Ar/l = -Ar/flfl (P 2 1) P 2 + \{v Q)q 2 P 6 

p — 1 p — I 2 \ 2 

+ i (p 2 - 30p + 90) q 3 (3 8 + -^{P 3 - 11 V + 1236 P - 252 °) 9 4 ^ 10 ^) A X 
+ ((fi 2 -l)(3 2 + qf3 6 (2m + p-6) + ^q 2 (3 s (l0m(p - 4) +p 2 - 30p + 90) 

+ iq 3 /3 10 (2m (l9p 2 - 276p + 630) + p 3 - lUp 2 + 1236p - 2520)^ A x 2 

1 — m ( 1 1 

+ — — f to/3 6 + - (p - 6)/3 6 + -g/3 8 (6m 2 + 6to(2 P - 9) + p 2 - 30p + 90) 

+ ig 2 /? 10 (30to 2 (3 P - 10) + 6?«(8p 2 - 129p + 308) + p 3 - 114p 2 + 1236p- 2520)) A x 

1 — 171 { 1 1 

+ m 2 /? 8 + - m (12p - 54)/3 8 + - (p 2 - 30p + 90) /3 8 

8 \ 6 6 

+q (^4m 3 f3 w + 4m 2 (bp - 19)/3 10 + i?«/3 8 (50p 2 /3 2 - 840p/3 2 + 2064/3 2 ) 
+ ^/3 8 (p 3 ^ 2 - H4p 2 ^ 2 + 1236p,3 2 - 2520/3 2 )^ A* 4 
' to 3 /3 10 + to 2 (5p - 19)/? 10 + im (50p 2 - 840p + 2064) /3 10 



10 V 24 



— (p 3 - 114p 2 + 1236p - 2520) /3 10 ^ A x 5 • (A2) 



I 

The full continuous free energy is defined via its ex- order at zero magnetic field read 
pansion coefficients in Eq. (|4T|) . Their values to the fifth 



.9(MH y, (A3) 
g(l,l) = 0, 1>1, (A4) 

ff (2,2) = i/3 4 (TO[0,0]-l) , (A5) 

«?(2,3) = -^ 4 to[0,1] , (A6) 

5 (3, 3) = ^ (/3 2 (to[0, 0] - 1)(2to[0, 0] +p - 6) 

+to[1,1]) , (A7) 
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g(2,4) = i/3 4 m[0,2] , (A8) 

g (3, 4) = ^/3 4 (2/? 2 m[0, 1] (4m[0, 0] + p - 8) 

+m(l,2)) , (A9) 

5 (4, 4) = i,/3 4 (2/3 4 (m[0, 0] - 1) (6p(2m[0, 0] - 5) 
96 v 

+6 (m[0, 0] 2 - 9m[0, 0] + 15) + p 2 ) 
+2/3 2 m[l, l](6m[0, 0] + p - 10) + m[2, 2]) , 



5 (2,5) = -/3 4 m[0,3] , 

3(3, 5) = ^/3 4 (3/3 2 (m[0, 2] (4m[0, 0] + p - 8) 



-4m [0, 1] 
1 



^1,3]) , 



(A10) 
(All) 

(A12) 



3(4, 5) = — /3 4 (6/3 4 m[0, 1] (6p(4m[0, 0] - 7) 

+6 (3m[0, 0] 2 - 20m[0, 0] + 24) + p 2 ) 

+3/3 2 (m[l,2](6m[0,0] +p- 10) 

+12m[0, l]m[l,l]) + m[2,3]) (A13) 



5(5,5) 



-B A (6^ 4 m[l, 1] (p(38m[0, 0] - 56) 



1440 r 

+36m[0, 0] 2 - 206m[0, 0] + p 2 + 212) 
+6/3 6 (m[0,0] - 1) (2p 2 (25m[0,0] - 57) 
+12p (I0m[0, O] 2 - 70?ti[0, 0] + 103) 
+24 (m[0, 0] 3 - 19m[0, O] 2 + 86m[0, 0] 
-105) + p 3 ) + 3/3 2 (m[2, 2] (8m[0, 0] + p 
-12) + 12m[l, l] 2 ) +m[3,3]) . (A14) 



The corresponding expansion coefficients of the order- 
parameter function to the third order are 



m[0, 0] 



p-2 



m[0,l]= -((3-2p)p + 6)/3 2 , 

m[l,l] = i(p(7p-24) + 12)/? 2 , 
m[0,2] = (p + 2)(p(p + 9)-27)/3 4 , 
m[l, 2] = -_(p(p(i7p + 19) - 228) + 204)/3 4 , (A19) 



m[2, 2] = -(p - 2)(p(25p - 32) - 44)/3 4 , 
m[0, 3] = -|(p(p((p - 29)p + 332) + 68) 



- 1504)^ 6 , 



(A15) 

(A16) 

(A17) 
(A18) 



m[l, 3] = -(p(p(2p(5p + 46) + 935) - 4296) 

+ 36 84)/3 6 , 
9 

m[2, 3] = --(p(p(p(52p + 285) - 1006) - 1068) 



(A20) 

(A21) 
(A22) 



3016)/3 6 
9 



(A23) 



m [3, 3] = -/3 b (p(p(p(283p - 488) - 1208) + 1184) 



-1776) 



(A24) 
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